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I .  INTRODUCTION 


Vol terra-series,  or  nonlinear  transfer-function  analysis, is  an 
efficient  and  highly  effective  method  for  analyzing  weakly  nonlinear  cir¬ 
cuits  under  small-signal  excitation.  It  has  been  used  to  analyze  nonlinear 
effects  in  a  wide  variety  of  circuits,  and  is  particularly  useful  for 
calculating  intermodulation  levels  in  MESFET  amplifiers. Volterra- 
series  techniques  have  also  been  used  to  analyze  time-varying  circuits  such 
as  mixers°^and  parametric  upconverters. '  However,  the  large  amount  of 
algebraic  effort  necessary  to  generate  analytical  forms  of  the  Volterra 
kernels  of  even  simple  circuits  has  prevented  the  wide  application  of 
Volterra-series  analysis.  This  analytical  complexity  has  forced  many  ex¬ 
perimenters  to  employ  simplifying  assumptions,  especially  in  the  topology 

H  o  t 

of  the  equivalent  circuits  of  solid-state  devices;  ’  these  approximations 
limit  the  accuracy  of  the  analysis.'- 

V 

This  report  describes  a  general-purpose  computer  program  that  performs 
a  Volterra-series  analysis  of  a  weakly  nonlinear  circuit.  The  program  is 
intended  primarily  for  use  in  the  design  of  microwave  circuits;  its  catalog 
of  circuit  elements  includes  the  distributed  elements  necessary  for  such 
work.  Because  the  program  formulates  and  solves  the  circuit  equations 
numericailly,  the  user  need  not  simplify  either  the  circuit  or  the  model  of 
the  solid-s'tate  device,  or  make  amy  of  the  other  cotmon  simplifying  assump¬ 
tions  often  found  in  Volterra-series  analyses  (e.g.,  that  the  input  fre¬ 
quencies  in  a  two-tone  excitation  are  nearly  identical). 
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II.  THE  NONLINEAR-CURRENT  METHOD 


The  program  uses  the  nonlinear-current  method and  a  nodal  formu¬ 
lation  of  the  circuit  equations.  In  this  approach  each  nonlinear  circuit 
element  is  described  as  a  linear  element  that  is  in  parallel  with  a  set  of 
current  sources;  each  current  source  represents  a  single  order  (greater 
than  1)  of  the  mixing  products,  and  its  current  is  a  nonlinear  function  of 
the  node-voltage  components  at  lower -order  mixing  frequencies.  The  weakly 
nonlinear  circuit  is  reduced  to  a  linear  circuit,  one  that  contains  the 
linear  elements  and  the  linear  parts  of  the  nonlinear  elements,  as  well  as 
a  set  of  excitation  sources. 

The  first-order  response  is  simply  the  response  of  the  linear  circuit. 
The  source  currents  at  second-order  mixing  frequencies  are  found  from  the 
first-order  node  voltages.  For  example,  in  the  case  of  a  simple  voltage- 
controlled  conductance  having  the  small-signal  I/V  characteristic, 

i  =  giv  g2V^  g3V^  ♦  •  ••  (1) 


where  the  first-order  terminal  ’oltage  v(t)  is 

,  Q 

=  2  I  V  exp(Ju»  t) 
q=-Q  ^ 

and  the  second-order  current^ is 


(2) 


l2(t)  =  g2Vi^(t) 


Q 

I 


q1  =  -Q 


J.,  'l.ql  h,q2 


(3) 


In  the  case  of  a  capacitor  having  the  charge /voltage  characteristic 


(4) 


=  C^v  +  C2V^  +  + 


the  second-order  current  i2(t)  is 


io(t)  = 


d[v^  (t)] 
2  dt 


-f  I  I 

^  q1=-Q  q2=-Q 


q1  "  ''l.qS  ®"PfJ^“q1  " 


One  can  find  the  second-order  node  voltages  by  solving  the  linear  network's 
admittance  equations  with  the  second-order  sources  used  as  excitations. 

The  process  is  repeated  to  find  the  third-  and  higher-order  node  voltages. 

Vol terra-series  analysis  is  considerably  more  efficient  than  gener¬ 
alized  harmonic-balance  analysis,  because  it  requires  neither  Fourier 
transformation  nor  iteration;  when  the  nonlinear-current  method  is  used,  • 

the  source  currents  are  found  directly  in  the  frequency  domain.  One  need 
not  evaluate  all  the  mixing  products  implied  by  Eqs.  (3)  or  (5);  it  is 
necessary  to  evaluate  only  the  small  set  of  products  that  contribute  to  the 
distortion  products  of  interest.  Thus,  if  only  the  minimum  number  of  sig¬ 
nificant  products  are  retained,  the  calculation  of  the  mixing  products  can 
be  performed  very  rapidly.  However,  Volterra-series  analysis  is  limited  to 
circuits  having  relatively  weak  nonlinearities;  circuits  having  strong 
nonlinearities  are  best  analyzed  by  means  of  harmonic-balance  techniques. 
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III.  PROGRAM  DESCRIPTION 


The  program  is  written  in  Turbo  Pascal  for  use  on  the  IBM  PC  or  PC- AT 
microcomputer.  The  program  is  menu  driven,  and  can  be  used  interactively 
to  generate  output  in  several  different  forms.  The  element  catalog  in¬ 
cludes  the  same  types  of  elements  found  in  most  linear -analysis  programs: 
linear  resistors,  capacitors,  inductors,  controlled  sources,  transmission 
lines,  stubs,  and  arbitrary  Impedances;  and  nonlinear  resistors,  capaci¬ 
tors,  and  controlled  sources.  At  present  the  program  is  limited  to  non¬ 
linear  elements  that  can  be  defined  as  voltage  controlled  (i.e.,  whose 
incremental  currents  are  single-valued  functions  of  voltage).  The  program 
has  a  modular  structure  that  allows  additional  types  of  circuit  elements  or 
additional  capabilities  to  be  added  easily. 

Figure  1  shows  the  program's  flow  chart.  After  reading  the  data  file, 
the  program  generates  the  frequency-mix  vectors^  of  only  the  desired  and 
contributory  mixing  frequencies;  i.e.,  it  ignores  mixing  components  of  the 
input  tones  that  do  not  contribute  to  the  distortion  product  of  interest. 
The  prograim  then  finds  the  source  currents  and  node  voltages  at  each  mixing 
frequency,  beginning  with  the  first-order  voltages  and  continuing  through 
the  order  of  the  mixing  product  of  interest.  At  each  order  the  program 
first  identifies  the  significant  lower-order  mixing  vectors  and  the  corre¬ 
sponding  voltages,  calculates  the  excitation  currents,  generates  the  nodal 
equations,  a."d  finally  calculates  the  node  voltages. 

Figures  2  through  A  illustrate  a  sample  calculation.  Figure  2  shows 
an  equivalent  circuit,  while  Figs.  3  and  4  show  the  data  file  and  output 
history,  respectively;  the  file  describes  a  single-stage,  narrowband, 
9.5-GHz  GaAs  MESFET  amplifier.  The  amplifier's  source  and  lead  impedances 
are  defined  by  tabular  data  in  the  files  ZSFET.DAT  and  ZLFET.DAT,  respec¬ 
tively  (the  program  interpolates  these  tables  if  they  do  not  contain  data 
at  precisely  the  correct  frequency).  The  FET  equivalent  circuit  includes 
the  three  nonlinear  elements  that  are  normally  dominant  in  establishing  the 
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•.  ttut  ******  t  **********  ************* 
**********  FET  IM  STUDY  ********** 

********************************** 

I 

FREQS  9.0E9  9.01E9 

dBm  0.0  0.0 

COEFS  -1  2 


CKT 

ZS 

1 

ZSFET.DAT 

ZL 

7 

ZLFET.DAT 

RES 

1  2 

2.0 

RES 

3  5 

2.0 

RES 

5  6 

1.5 

RES 

4  7 

1.5 

IND 

6  0 

0.04E-9 

CAP 

2  4 

0.030E-12 

CAP 

4  5 

0. lE-12 

VCN 

2  3 

4  5  0.03311 

0.00401 

-0.00426  0  0 

DIDDEC 

2  3  0.45e-12 

0.7 

0.500  -1.5 

GNL 

4 

5  0.00367  - 

0.000370 

0.000144  0  0 

END 


Fig.  3.  Data  File 
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Pile  Nate:  :eti«.aat 


laitiai  Inpat  Freos  and  peters 
Coiponent  Freo.  PoteridBii 


*  1 

Q.;'E*000'? 

O.OOf't) 

i  ■' 

a.C-E^OOO^ 

0.0000 

pf  Id  praduct: 

MN  =  2 

(2  -  71 

INPUT  FRE8UENCIES. 

6H: 

n 

f2 

f3 

14 

«. 000000 

4.010000 

0.00 

0.00 

?. 100000 

4.110000 

0.00 

0.00 

N.IOOOOO 

4.210000 

0.00 

0.00 

-.300000 

4.310000 

7.00 

0.00 

=.400000 

4.410000 

0.00 

0.00 

4.500000 

4.310000 

0.00 

0.00 

4.600000 

4.61000':^ 

0.00 

0.00 

4.700000 

4.710000 

0.00 

0.00 

4.800000 

4.810000 

0.00 

0.00 

4.400000 

4.410000 

0.00 

0.00 

10.000000 

10.010000 

0.00 

0.00 

Order  of  !tl  prpduct  =  ’ 

FI  OUTPUT  PARAI1S 


<5 

ReUs) 

UlZs) 

IiiZD 

0.00 

10.30 

45.60 

100.00 

6.00 

0.00 

10.30 

45.60 

100.00 

6.00 

0.00 

10.30 

45.60 

100.00 

6.00 

0.00 

10.30 

45.60 

100.00 

6.00 

0.00 

10.30 

45.60 

100.00 

6.00 

0.00 

10.30 

45.60 

100.00 

6.00 

0.00 

10.  JO 

45.60 

100.00 

6.00 

0.00 

10.30 

45.60 

100.00 

6.00 

0.00 

10.30 

45.60 

100.00 

6.00 

0.00 

10.30 

45.60 

100.00 

6.00 

0.00 

10.30 

45.60 

100.00 

6.00 

Pout 

in  FREQ.  6H: 
iin 

I  in  OUTPUT  PRRUnS 
Re(ZL)  li(ZL)  Pout 

10.27 

4.020000 

100.00 

6.00 

-17.53 

10.23 

4.120000 

100.00 

6.00 

-17.70 

10.18 

4.220000 

100.00 

6.00 

-17.88 

10.13 

4.320000 

100.00 

6.00 

-18.06 

10.07 

4.420000 

100.00 

6.00 

-18.26 

10,01 

4.520000 

100.00 

6.00 

-18.46 

4.44 

4.620000 

100.00 

6.00 

-18.67 

4.87 

4.720000 

100.00 

6.00 

-18.88 

4,80 

4,820000 

100.00 

6.00 

-14.10 

4.72 

4.420000 

100.00 

6.00 

-14.32 

4.64 

10.020000 

100.00 

6.00 

-14.54 

Fig.  4.  Output  Listing 
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device's  intermodulation  (IM)  characteristics:  the  gate/source  capaci¬ 
tance,  the  transconductance,  and  the  drain/source  resistance. 

The  data  file  lists  the  excitation  frequencies,  the  power  levels,  and 
the  coefficients  of  the  IM  product  of  Interest,  followed  by  the  circuit 
description.  The  elements  VCN,  GNL,  and  OIODEC  are  nonlinear;  VCN  and  GNL, 
the  controlled  current  source  and  drain/source  conductance,  are  described 
by  the  power-series  coefficients  of  their  incremental  I/V  characteristics. 
DIODEC,  representing  the  gate/source  capacitance,  is  the  capacitance  of  an 
ideal,  uniformly  doped  p-n  Junction.  The  capacitance  is  described  by  its 
Junction  parameters  (Cjq,  ♦,  and  y)  and  bias  voltage;  the  program  calcu¬ 
lates  the  coefficients  of  its  charge/voltage  characteristic  from  these  data. 

Figure  4  shows  the  results  of  a  two-tone  sweep  across  the  amplifier's 
passband:  it  lists  output  powers  of  the  linear  responses  at  the  excitation 
frequencies  f^  and  fg,  as  well  as  third-order  IM  products  at  the  mixing 
frequency  2f2  -  f ^ .  The  output  listing  includes  the  source  and  load 
Impedances  at  the  input  and  IM  frequencies.  Computation  speed  is  limited 
primarily  by  the  solution  of  the  admittance  equations,  and  is  almost  inde¬ 
pendent  of  the  number  of  nonlinear  elements;  thus,  many  more  elements  could 
be  treated  as  nonlinear  ones  without  a  significant  increase  in  calculation 
time.  The  data  in  Fig.  4  require  only  2.5  seconds  of  computer  time  per 
line. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect -engineer"  for 
national  security  projects,  specializing  in  advanced  military  space  systems. 
Providing  research  support,  the  corporation’s  Laboratory  Operations  conducts 
experimental  and  theoretical  invest igat ions  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  Is  the  technical  staff’s  wide-ranging  expertise  and  Its 
ability  to  stay  current  with  new  developments.  This  expertise  Is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 

Aerophysics  Laboratory:  Launch  vehicle  and  reentry  tluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  excimer  laser  development  including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out-of-f  leld-of -view  rejection,^ 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atonic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  Intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications; 
microwave  semiconductor  devices,  microwave/millimeter  wave  measurements, 
diagnostics  and  radlometry,  raicrowave/milllmeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  In  space  and  enemy-induced 
environments. 

Space  Sciences  Laboratory:  Magnetospherlc,  auroral  and  cosmic  ray 
physics,  wave-parClcle  interactions,  magnetospherlc  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomy, 
Infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation. 


